On Conserved Current in Markovian Open Quantum Systems 
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We reexamine the markovian approximation of local current in open quantum systems, discussed 
recently by Gebauer and Car. Our derivation is more transparent, the proof of current conservation 
becomes explicit and easy. 



I. INTRODUCTION 

Quantum mechanics is reversible at the microscopic 
level. If, however, an atomic system is coupled to a 
reservoir the atomic reduced dynamics becomes irre- 
versible and the corresponding evolution equation con- 
tains a memory term. Yet, in Markov approximation the 
time evolution of the atomic reduced density matrix p 
satisfies the memory-less master equation: 



Cp, 



(1) 



where H is the atomic Hamiltonian. The irreversible 
term can be written into the Lindblad-form 0,@: 



(2) 



where {V^} are certain atomic operators and {<?;} are 
complex coefficients satisfying Re gi > 0. This equation 
constitutes the phenomenological dynamics of the given 
markovian open quantum system. 

In Hamiltonian dynamics with local potentials and in- 
teractions an important conserved quantity is the particle 
density. Consider, for simplicity's sake, a single atomic 
electron with the canonical pair of operators x, p. The 
electron's density operator at location r reads: 



n(r) = S(r — x) = \r) (r\ 



(3) 



where the position eigenstates satisfy the equation 
x |r) = r \r). There is a local current defined by: 



J(r) 



1 

2m 



\pn{r) + h(r)p] 



(4) 



The current j is conserved, i.e., the Heisenberg operators 
fit, J t satisfy the continuity equation: 



dh t (r) 
dt 



VJ t (r) = 0. 



(5) 
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A similar continuity equation follows for the expectation 
values as well. The continuity equation still holds if the 
atom interacts with other systems provided the inter- 
action is local. Let us, indeed, suppose a certain local 
interaction with a reservoir. The atom becomes an open 
system governed by the master equation Q in the marko- 
vian approximation. Let us write down the continuity 
equation for the expectation values of n and J in the 
atomic state p, and substitute eq.Q. In Schrodinger 
picture we obtain the following: 



d(n(r)) 
dt 



V(J(r)> = (C*n(r)). 



(6) 



We see that the extra term (C*n(r)) = tr[n(r)Cp] on the 
r.h.s. may in general violate the local conservation of the 
current. Gebauer and Car [3| noticed that such violation 
is a consequence of Markov approximation. They also 
pointed out that current conservation can be restored by 
adding a dissipative correction J d to the Hamiltonian 
current J. These authors have derived J d from the ex- 
act reversible dynamics of the atom+reservoir. Here we 
re-consider the issue and propose a shorter and trans- 
parent derivation. Our formulating Jd is fairly explicit 
to satisfy local conservation. This latter is the main 
progress with respect to ref.Q which would need sub- 
tle and lengthier proof 0, not even published in ref.Q 
or elsewhere. 

It is clear from eq.© that J + J d will satisfy the 
continuity equation if: 



VJ D {r) = C*n(r). 



(7) 



In one dimension with boundary conditions 
«7d(±oo) = the above equation has a unique so- 
lution J o{ r ) = /^oo dr 1 C*h{r') but in higher dimensions 
additional considerations are necessary to uniquely 
determine J d- In fact, the Markov dynamics (JIJ of the 
atomic state p is not sufficient to calculate J d- We must 
carefully inspect the Markov approximation of the exact 
current. In Sec. II we review the derivation of the master 
equation Q in the usual Born-Markov approximation. 
Then Sec. Ill contains our derivation of the dissipative 
current J rj. 



2 



II. BORN-MARKOV APPROXIMATION 

Let stand H+K+Hr for the total Hamiltonian where 
H belongs to the atomic system, H R belongs to the reser- 
voir. The interaction Hamiltonian is local, we consider 
the simple case K = V(x)F where the "field" F is a 
reservoir operator 0. We suppose that the initial state 
of the whole system is pp R where p R is the reservoir equi- 
librium state. We denote by O t the operator O at time 
t in interaction picture. Let us introduce the reservoir 
correlation function 



g{a - t) = tr R (^F a F T p R 



(8) 



which is a complex non-negative time-translation invari- 
ant kernel. Consider the change of the atomic reduced 
density matrix during a certain time A. In interaction 
picture, second order perturbation theory yields: 



PA 



1 



A 



da ( dr (v T pV a - VaV T pj g(a - r) 



JO 



+H.C. (9) 



where the first order term in V vanishes since we assume 
tr R (Fp R ) = 0. Suppose that ta S> t r where ta is the 
time scale of atomic state evolution and t r is the reservoir 
correlation time. Choose A in between the reservoir and 
atomic time-scales: 



t r < A < ta- 



(10) 



A standard way to proceed is the Born-Markov approx- 
imation that requires weak coupling and results in the 
master equation Q for the atomic state pt time-coarse- 
grained on scale A. 

The implementation of the Born-Markov approxima- 
tion has two steps. We identify the l.h.s. of by the 
time derivative dp/dt of the time-coarse-grained atomic 
state p. On the r.h.s., we extend the integration limit as 
well as the denominator A to oo. The latter limit is rigor- 
ously justified by a particular rescaling [6j where V — > 
while V 2 A is constant 0- Accordingly, the explicit form 
of the Lindblad-superoperator becomes: 



Cp 



lim — 

t^oo t 



t 

da 

o 



dr 
o 



(VrpV a 



V a V T p)g(a-T) 

+H.C. (11) 



It is easy to inspect that the r.h.s. is already of the 
Lindblad form (J2J. To achieve the familiar spectral rep- 
resentation, we expand the interaction potentials as 



(12) 



where the sum is over all possible transition frequencies 
w = e n — £m of the system Hamiltonian H = ^ n e n \n) (n|. 
Obviously, the spectral component of V is of this form: 



E 



\n)V n 



(13) 



As t — > oo the oscillating terms in the integral l|ll|) drop 
out and we get the master equation where: 

£p = E(^^ - v^p)gi +) + hc- (i4) 



and 



,(+) 



dre^g(r). 



(15) 



III. 



THE MARKOVIAN CONSERVED 
CURRENT 



Constructing the correct markovian coarse-grained 
current J + Ju needs special care. First, we calculate the 
average of the exact current over the period A of time- 
coarse-graining. We expect that it becomes the sum of 
the unperturbed contribution J and the markovian cor- 
rection J f) in the Born-Markov limit: 



1 



A 



da{J) a 



(J) + (Jd) 



(16) 



where (. . . ) (T stands for expectation value in the ex- 
act time-dependent state evolved from pp R by the to- 
tal Hamiltonian H + H R + K while (...) stands for the 
expectation value in the coarse-grained state p. Let us 
introduce the time-dependent auxiliary variable Z(r;t) 
defined as 



Z(r;t) = / daJ a {r) 



(17) 



Then, in the interaction picture, the l.h.s. of l|ltjfl can be 
written as: 

^trZ(A)p A - ltr^ A daZ(a)^. (18) 

In Born-Markov limit, the time-coarse-grained atomic 
state pt evolves smoothly and only "differentially" during 
time A, so does the current J a , too, in such states: we 
can thus write the first term as just (J). In the second 
term, before taking the Born-Markov limit, we substitute 
the perturbative expression of dp a j da: 



1 f A - 
-— tr / daZ(a) 
A Jo ■'<> 

-H.C. 



dr 



(VrpVa 



V a V T p)g(a-T) 
(19) 



Now we take the Born-Markov limit which, again, al- 
lows us to calculate the r.h.s. by taking the limit A — » oo: 

(Jd) = 

- lim -tr [da fdrZ(a) (v T p% - V„V T p\ g(a - r) 

t^oo t J J \ J 



H.C. 



(20) 
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We can express the dissipative current itself: 
Jd = 

- lim - (do [dr (v a Z(a)V T - Z(a)V a V T ) g{a 

S->qo t J Q J \ J 

+ H.C. 



(21) 



Recall that Z(r;a) = f° d\J\(r) so that J d{t) has 
turned out to be a local linear functional of the unper- 
turbed local current J{r). 

It is easy to prove that J d contributes properly to 
current conservation. Let us take the divergence of both 
sides of eq. (|2*TT) . Observe that 



VZ(o-) = n a — h, 



(22) 



because density conservation VJ ff + dfi a /da = holds 
for the unperturbed atomic current and density. We have 
assumed that the interaction is local, so [rio-jV^] = 0. 
Therefore the terms proportional to fi a cancel each other 
on the r.h.s. of ea. (|21|l and we obtain: 



VJ 



(23) 



t^oc t 



lim - da dr V a hV T - nV a V T g(a - r) + H.C. 



The r.h.s. coincides with C*n as it should, cf. eq.0). 



IV. CLOSING REMARKS 

Our central result is the expression l|21l) of the dissi- 
pative current Jd- Its form was convenient when we 
proved the continuity equation from the locality of the 
interaction. In concrete calculations, however, the spec- 
tral representation is the preferred one. As the spectral 



representation of Z(a) is 



iujj 



(24) 



ea. (|21|l yields two terms for the spectral expression of 
J d- To compare it with the result of Car and Gebauer 
[3( , we present the spectral expression of the expectation 
value: 



/ 



(Jd) = tr 



E 



E 



(25) 



U)J-\-UJ\— LJ2— LJ 1 — UJ2 — . 

ui p —0 ujj-\-ujp=0 J 



Here J Uj ,Pu p and V u are the respective spectral repre- 
sentations of J,/5 and V, cf. ea. l(T^|) . The two sums 
cancel the singularity at uij = 0. We note without go- 
ing into the details that, for non-degenerate transition 
frequencies, our result l|25|l coincides with eq.(16) in Q 
apart from minor typos of the latter Q . 

Our result allows us to calculate the markovian cor- 
rection Jd to the local current in open quantum systems 
with discrete spectrum. The method applies to quantum 
dots directly [|| ■ The case of quantum brownian motion, 
however, requires a suitably modified approach to cope 
with the markovian limit of an open system with contin- 
uous spectrum. 
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